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ABSTRACT 

Clumps in the solar neighborhood's stellar velocity distribution could be caused by spiral density 
waves. In the solar neighborhood, stellar velocities corresponding to orbits that are nearly closed 
in the frame rotating with a spiral pattern represent likely regions for stellar concentrations. Via 
particle integration, we show that orbits can intersect the solar neighborhood when they are excited 
by Lindblad resonances with a spiral pattern. We find that a two-armed spiral density wave with 
pattern speed placing the Sun near the 4:1 Inner Lindblad Resonance (ILR) can cause two families of 
nearly closed orbits in the solar neighborhood. One family corresponds to square shaped orbits aligned 
so their peaks lie on top of, and support, the two dominant stellar arms. The second family correspond 
to orbits 45° out of phase with the other family. Such a spiral density pattern could account for two 
major clumps in the solar neighborhood's velocity distribution. The Pleiades/Hyades moving group 
corresponds to the first family of orbits and the Coma Berenices moving group corresponds to the 
second family. This model requires a spiral pattern speed of approximately 0.66 ± 0.03 times the 
angular rotation rate of the Sun or 18.1 ± 0.8 km s _1 kpc -1 . 

Subject headings: Galaxy: kinematics and dynamics — Galaxy: disk — stars: kinematics 



1. INTRODUCTION 

The velocity distribution of stars in the solar neigh- 
borhood contains structure which has been particularly 
clearly revealed from recent studies of Hipparcos obser- 
vations jDehnen|[rp98: Skuli a^et^LJIl^Q^jF^ma^ve^^J 
l27^IChereul et alJll998L ll999HNordstrom et alJl2004(1 . 
Much of this s tructure was previously identified with 
moving groups (Eggen 1996). Moving groups are groups 
of stars which are kinematically associated. Young, 
early-type stars can be moving together becaus e they 
carry the kinematic signature of their birth (Eggcn 
1996). However, a number of the kinematic clumps 
identified through the kin ematic studies also contain 
later t ype and older stars llDehnenl 119981 IFamaev et al.l 
120041 iN ordstro m et al.l I2004D. ^ Dubbed superclusters 
ijFamaev et all 120041 lEggenl 11996(1 . these kinematic as- 
sociations of stars, which are observed all over the sky, 
can share the same space motions as well-known open 
clusters. The best documented of the superclusters are 
the Hyades supercluster associated with the Hyades clus- 
ter, the Sirius or Ursa Major supercluster, associated 
with the Ursa Major cluster, and the Pleiades moving 
group (or Local Association) associated with young clus- 
ters such as the Pleiades. 

There have been a number of studies ex ploring the ori- 
gin of moving groups and superclusters. lEggenl 1(19961) 
and others have proposed that the structure in the ve- 
locity distribution was a result of inhomogeneous star 
formation and the dissolution of clumps of stars formed 
simultaneously. Large moving groups could be produced 
by the dissolution of large stellar agglomerat ions associ- 
ated with spiral arms (e.g.. lAsiain et al.lll999(l . However, 
the presence of older stars in the moving groups or su- 
perclusters has presented a challenge for this simplest 
scenar io. To account for the older stars iChereul et al.l 
( 1998) suggested that superclusters could be superposi- 



tions of different age moving groups. Because some of 
the moving groups contain stars spanning a range colors, 
groups could be long lived and may not necessarily have 
transient features associated wi th recently formed clus- 
ters that are being disrupted. IDe Simone et alJ (2004) 
suggested that the clumps in the velocity distribution are 
due to irregularities in the Galactic potential. In their 
model, spiral arms near corotation were induced stochas- 
tically on a distribution of stellar test particles causing 
localized structure in the stellar velocity distribution. 

The possibility that we consider here is that the clumps 
in phase space (in the velocity distribution) are associ- 
ated with or caused by spiral density wav es traveling in 
the so l ar neighborhoo d, as discussed by IFamaev et alJ 
l(2004|) : IDehnenI l(l99^ : IDe Simone et all l(2004() . This 
type of explanation would not be unprecedented. The 
Hercules stream at a tangential velocity of v ~ 45km/s 
has been explain ed in terms of perturbat i ons caused by 
the G alactic bar l(Dehnedl200TilFuxll200lHRaboud et alJ 
1998). This stream is due to stars in orbits strongly 
affected by the Galactic Bar's 2:1 Outer Lindblad Reso- 
nance. Based on this elegant explanation for the Hercules 
stream, we are motivated to explore similar explanations 
for structure in the velocity distribution, but at velocities 
nearer those of circular orbits. Before we consider the ef- 
fect of spiral density waves on the solar neighborhood's 
velocity distribution, we first review what is known about 
the strength, number of arms and possible pattern speed 
(angular rotation rate) for spiral density waves likely to 
be propagating in the s olar nei g hborh ood. 

The recent study of iValled ((2002) provides a good 
summary of the many studies that have used observa- 
tions to map the Milky Way disk. Cepheid, HI, CO 
and far-infrared observations suggest that the Milky Way 
disk co ntains a four armed tig htly wound structure, 
whereas (Drimmel fc Spergel 1(200 1() have shown that the 
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near-infrared observations are consistent with a domi- 
nant two-armed structure. The dominant two-armed and 
weaker four-armed struc ture was previously proposed by 
lAmaral fc Lepinel l)1997j) . The nearest spiral arm (ex- 
cluding the local Orion armlet) is the Sagittarius/Carina 
arm, 0.9kpc away from the Sun in the direction toward 
the Galactic center. The distance between this arm and 
the Perseus arm in the opposite direction from us (to- 
ward the Galactic anticenter) is about 2.5kpc. 

Based on observational constraints on the Galactic spi- 
ral structure, a number of studies have created dynam- 
ical models to fit kinematic observations. These mod- 
els are sensitive to, and so constrai n, the spiral p attern 
speed. Reviewing previous work, iShavivl ((2003) finds 
a clustering of estimates for the pattern speed of lo- 
cal spiral structure near Q s ~ 20kms _1 kpc , though 
other studies suggest £l s ~ 13kms _1 kpc . iLenine et alJ 
( 2001) suggest that locally the Milky Way can be mod- 
eled by the superposition of a two- and four-armed struc- 
ture. Their model places the Sun near the corota- 

tion resonance fl s ~ 28kms _1 kpc ), and was fit to 
Cep heid kinematics. The recent gas dyna mical stud- 
ies l)Martos et alJl2005t iBissantz etafll2003|) match the 
properties of the gas in nearby arms with a spiral pat- 
tern speed of ~ 20kms _1 kpc . iMartos "eTal] (2005) pro- 
pose that a two-armed stellar structure consistent with 
the stellar distribution inferred from COBE could cause 
four-arms in the gas distribut ion near the Sun. Th e gas 
dynamical model proposed by Bissantz et al. (2003) with 
a similar spiral pattern speed matches HI and CO kine- 
matics. 

In this paper we consider the effect of a spiral density 
wave on the solar neighborhood velocity distribution. In 
section 2, we discuss how the velocity components of stars 
near the Sun are related to their orbits in the Galaxy. In 
section 3 we describe our technique for determining the 
velocity components of different populations of stars. In 
section 4 we explore how different spiral density waves 
perturb the velocity distribution. A summary and dis- 
cussion follows. 

2. EPICYCLIC MOTION AND THE POSITION ON THE U, V 

PLANE 

We must relate the observed velocities of stars to quan- 
tities describing their orbital motion in the Galaxy. The 
solar neighborhood velocity distribution can be described 
as a function of the azimuthal velocity component, v, and 
the radial velocity, v r = —u, where u > corresponds to 
velocities toward the Galactic center. We define v such 
that the tangential component of the velocity, in the di- 
rection of Galactic rotation, is equal to Vq + v. where 
Vq = 220km/s, is the local standard of rest at the posi- 
tion of the Sun fo r a Galactocentric radius of Ro — 8kpc 
ijReid et 3.1.1(19991) . 

In the absence of perturbations from spiral arms, the 
motion of stars in the disk of a galaxy can be described 
in terms of radial or epicyclic oscillations about a cir- 
cular orbit. It is useful to specify the relation between 
the observed velocity components u, v and parameters 
that describe the epicyclic motion. These parameter are 
the mean radius or guiding radius r g and the epicyclic 
amplitude. The energy of an orbit in the plane of an ax- 
isymmetric system (neglecting perturbations from spiral 



structure) is 
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where the potential energy, lnr, is that appropriate for 
a flat rotation curve, and r is the Galactocentric radius. 
In the above equation we have put velocities in units of 
Vo and radii in units of i?o = 8.0 kpc, the distance of the 
Sun from the Galactic Center. The angular rotation rate 
of a star in a circular orbit at the Sun's radius from the 
Galactic Center, Qq = 28km s _1 kpc -1 which is based on 
obser vations of the proper motion of Sag A* l)Reid et alJ 
1999) and is consistent with me asurements based on Ty - 
cho and Hipparcos observations ( Oiling & Dehncn 2003) . 
Henceforth we place angular rotation rates and pattern 
speeds in units of £Iq. 
In an epicyclic approximation we can write the energy 
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where \ + In r g is the energy of a star in a circular orbit 
about a guiding radius r g and E ep i is the energy from 
the epicyclic motion. 
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where a is the epicyclic amplitude and n is the epicyclic 
frequency at the guiding or mean radius r g . We can also 
write E epi — jn where j is the radial action variable. 

We now consider stars specifically in the solar neigh- 
borhood restricting us to a specific location in the 
Galaxy. The velocity distribution in the solar neighbor- 
hood, or number density of stars as a function of u, v is 
measured near the Sun where r ~ 1 in units of Rq. Us- 
ing r = 1 in equation Q and setting the energy equal to 
that written in terms of the epicyclic motion using equa- 
tions EJ , we solve for r g . It is convenient to define 
the distance between the guiding or mean radius and the 
Sun's Galactocentric radius, s = r g — 1. To second or- 
der in v we find that a star near the Sun with velocity 
components u, v has a guiding radius with 

s « v (4) 



and epicyclic amplitude 



(5) 



The above two relations allow us to relate the veloc- 
ity components u, v for stars in the solar neighborhood, 
to quantities used to describe the epicyclic motion; the 
guiding radius and epicyclic amplitude. Particles with 
positive v have s > and so guiding radii that are larger 
than Rq, and mean radii outside the Sun's Galactocen- 
tric radius. Particles with negative v have s < and 
so are expected to spend most of their orbits inside Rq. 
The distance from the origin u = v = determines the 
epicyclic amplitude. 

The angular rotation rate of a star is most sensitive to 
the value of its guiding radius r g . The location of Lind- 
blad resonances with a periodic perturbation are there- 
fore primarily set by a star's guiding radius. In the solar 
neighborhood the guiding radius is primarily dependent 
on the v velocity component. The location of resonances 
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on the u, v plane is therefore most sensitive to the v value. 
We focus on the location of resonances because they are 
one possible cause for the structure in the solar neighbor- 
hood velocity distribution. In the vicinity of a resonance 
there is a bifurcation in the families of periodic orbits and 
there are no orbits that have low epicyclic motion (e.g., 
lContopoulosll975j) . Because resonances can induce large 
epicyclic (or radial) motions, stars with guiding or mean 
radii distant from the solar neighborhood could be seen 
in our vicinity, as is the case for the Hercules stream and 
the Galactic Bar's Outer Lindblad Resonance. 

3. FINDING ORBITS THAT ARE NEAR PERIODIC ORBITS 

In this section we describe our numerical procedure 
for constructing a function that can be compared to the 
number density of stars f(u,v)dudv in the solar neigh- 
borhood. We are restricting this initial study to a spiral 
perturbation moving at a unique fixed pattern speed or 
angular rotation rate f2 s . In this case the Jacobi integral 
is conserved. By transferring to the frame rotating with 
the pattern, the Hamil tonian can be brought to a time 
independent form (e.g.. lBinnev fc Tremainelll987[) . 

We wish to weight positions in the u, v plane at the 
position of the Sun according to how likely they are to 
be occupied by stars. Since stars are born in spiral arms 
with low velocity dispersion, we assume that they are 
born in orbits that are nearly periodic or closed in the 
frame rotating with the arms. After a star is born, it 
can be scattered by molecular clouds or transient spiral 
structure as it orbits in the Galaxy. Older stars would 
occupy orbits that oscillate at a greater extent about the 
periodic or closed orbits. We have chosen to weight or- 
bits according to how near they are to closed or periodic 
orbits. Consequently we construct a numerical way to 
measure the extent of oscillation about a periodic orbit 
for stars in the solar neighborhood. 

For each position in the u, v plane, we carry out a nu- 
merical integration of a particle with initial condition 
r = 1,0 = corresponding to its location near the Sun. 
Its initial radial velocity v r = — u and tangential veloc- 
ity 1 + v (in the inertial frame). In the frame rotating 
with the spiral pattern, we integrate the trajectory of the 
particle. Each time the particle passes through — 
(the Sun-Galactic center line), we record its position and 
velocity components. After a specified number of or- 
bits (in most cases 10), we compute the variance of radii 
(ay) and u values (ay) at times when the particle passed 
through = 0. The weighting function we used to esti- 
mate the likelihood of populating the position in the u, v 
plane is the sum of these two variances or W = ay + ay . 
We adjusted the number of orbits integrated so that the 
resulting weighting function was fairly stable; doubling 
the number of orbits integrated did not result in large 
changes in W. By keeping the number of integrated or- 
bits low, we could explore a larger region of parameter 
space. 

We use our weighting function to identify locations on 
the u, v plane likely to be populated with stars. Orbits 
with initial u, v values that have low values of our weight- 
ing function W correspond to orbits where a young star 
is likely to be born. Older stars could remain in these 
locations but also would also be scattered out to larger 
values of W corresponding to locations in the u, v plane 
with fewer stars. 



Our procedure for constructing a weighting function 
is related to that used to construct surfaces of section. 
Surfaces of section can be computed by fixing the Ja- 
cobi integral and plotting u vs r each time the par ticle 
passes through = fe.g.. lBinnev fc Tremainell987D . In 
a surface of section fixed points correspond to periodic 
or closed orbits in the galaxy. In such a diagram, our 
weighting function describes the square of the distance 
in phase space from a fixed point. We do not use surfaces 
of section here because they are generated for each value 
of the Jacobi integral. We wish to study the likelihood 
of populating orbits as a function of velocity in the solar 
neighborhood so we desire a weighting function that is a 
function of u, v. 

Our method can be co mpared to the b ackwards inte- 
gration technique used bv iDehnenl {2000). Stars at each 
u, v were integrated backward in time while the Galactic 
bar was reduced in strength. Velocities on the u, v plane 
corresponding to orbits with initially low epicyclic mo- 
tions (b efore the grow th of the Bar) were given higher 
weights (|Dehnen| [20001. The theory of adiabatic invari- 
ants implies that particles with initially low epicyclic mo- 
tion settle onto periodic or closed orbits following the 
slow growth of a perturbation. Consequently, we ex- 
pect that the backwards integration method produces 
weighting contours similar to those of our weighting func- 
tion. Compared to the backwards integration method 
our weighting technique has the advantage that it is in- 
sensitive to the initial stellar velocity distribution and the 
manner of perturbation growth (see iFuxl 120011) . These 
assumptions are a particular problem for spiral density 
waves si nce they could vary in b oth amplitude (as consid- 
ered bv lDe Simone et apl2004[) and wa vevector or pitch 
angle (as considered bv lFuchsll2001alHL 

Our particle integration is done for particles moving in 
the plane of the Galaxy. We assume a gravitational po- 
tential V(r, 0, t) = Vo(r)+Vi(r, 0, t) and Vb(r) = lnr cor- 
responding to a flat rotation curve. We assume that the 
gravitational potential perturbation Vi(r, 0, t) is caused 
by tightly wound logarithmic spiral density perturba- 
tions. Here r and are the radius and azimuthal angle 
in the Galactic plane and t is time. The potential pertur- 
bation can be expanded in terms of Fourier components 

Vy(r,0,i) = J] A m cos(a ro In r + m(0 - Q a t - 7m)) ■ (6) 



If only one Fourier component is present, m is the number 
of arms of that component. The angle 7 m corresponds 
to an angular offset measured at time t — 0. Measured 
at the Galactic Center, 7 m is the angle between the lo- 
cation of the Sun (r = 1,0 = 0) and the peak of the 
spiral pattern at the same radius but at an azimuthal 
angle differing from that of the Sun. A subsequent figure 
illustrates this angle. The parameter a m depends on the 
pitch angle, p m where a m — mcotp m . Note a m < for 
trailing arms when the rotation is clockwise. Maps of the 
Milky Way are commonly shown from the view point of 
an observer located above the Galaxy in the direction of 
the North Galactic Pole (e.g., iDrimmel fc Spergell2001l: 
lValleejl2002l On these maps the rotation is clockwise. 

In the WKB or tight winding approximation, the am- 
plitude of the potential perturbation Fourier component 
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is related to the density perturbations 

— 2TtGEoSmRo 
m ^ i iT72 v) 

KM K 

(e.g.. iBinnev fc Tremainelll987|) . In the above equation 
we have placed A m in units of Vq. Here E m is the 
amplitude of the to— th Fourier component of the spi- 
ral mass surface density and S m = E m /£o- The mean 
surface density of d isk mass in the solar nei ghborhood 
is £ ~ 50M Q pc- 2 l|Holmberg fc FlvniJ l2000). We esti- 
mate for the solar neighborhood 

S \ /220kms -1 ^ 2 ( 7 \ f R Q 



A 



0.035. 



V 50M©pc- 



V 2 



V8kpc 

A two-armed spiral is seen in the near infrared 
COBE/DIRBE data with pitch angle in the range p ~ 
15.5 — 19°, corres ponding to a = to cot p i n the range 
5.8-7.2 for to = 2 iforimmel fc Snergelll2fioT|) . 

We now estimate ballpark values for the potential 
perturbation component A-x- Usin g 52 ~ 0.15 consis- 
tent w ith parameters described bv IDrimmel fe Spergell 
(|2001l) for the stellar component (based on K band 
observations) the above equation |H1 gives us A 2 

lergell 



0.005. Specifically Dri mmel fc Spergell l)2001|) estimated 
(Smaa — Smin)/S m i7, = 0.32, though they suggested that 
the true value could be higher. 

4. STRUCTURE IN THE U, V PLANE 

Because the Milky Way has a dominant two-armed spi- 
ral pattern, we consider dominant two-armed spiral den- 
sity waves. In Figure ^ we show the structure of our 
weighting function W(u, v) for a pure cosine two-armed 
spiral pattern with parameters listed in Table ^ Here- 
after Tabled lists parameters used to describe the spiral 
perturbation. Each line in the table is labeled with a cor- 
responding figure displaying the output of a simulation. 

The strength of the spiral pattern of the simulation 
shown in Fig ^ is near but somewh at above that esti- 
mated bv IDrimmel fc Spergel l)2001() with A 2 = -0.005. 
Peaks in density correspond to dips in the gravitational 
potential. We use negative Fourier amplitudes so that 7 m 
refers to the location of the peaks (with respect to the 
Sun) of the spiral density pattern. The density perturba- 
tion has the opposite sign of the potential perturbation. 

Each panel in Fig ^ corresponds to spiral structure 
moving at a different pattern speed. We now relate po- 
sitions in the u, v plane to those expected for Lindblad 
resonances. The to : 1 Inner Lindblad resonance is lo- 
cated where _ 

Q s = Q(l-V2/m). (9) 

for a flat rotation curve, where k = ^/2fl and ft and k 
are computed as a function of the guiding radius. Here 
we adopt the notation to : 1 ILR referring to to epicyclic 
oscillations for each orbit around the galaxy in the frame 
moving with the pattern. Specifically, for the 4:1 ILR 
n s = 0.65ft, for the 3:1 ILR ft s = 0.53ft, and for the 2:1 
ILR ft s = 0.29ft. The topmost panel of Figure □ with 
ft s = 0.4 (in units of fto) has the Sun outside but near 
the 2:1 ILR. The bottom panel (with ft s = 0.7) has the 
Sun just outside the 4:1 ILR and so nearer to corotation. 

In the images shown in Figure ^ the log of our weight- 
ing function, log 10 W(u, v) is shown in gray scale. Con- 
tours are shown with the highest contour with log 10 W = 



— 1 and a contour interval of 1. The blackest regions 
correspond to smallest values of W with log 10 W ~ — 4 
which are those most likely to be populated by young 
stars. These regions could also maintain older stars at 
low velocity dispersion. 

In Figure ^ structure is primarily influenced by the 
value of w, as expected from our discussion in section 
2; since v sets the location of resonances with the spiral 
pattern. For larger values of v the angular rotation rate 
about the Galaxy drops. As v varies, resonances with the 
spiral perturbation are crossed. Consequently orbits af- 
fected by resonances that have large epicyclic excursions 
.can be seen in the Solar neighborhood. 

Dark regions on these diagrams correspond to nearly 
closed orbits in the frame moving with the spiral pattern. 
The strongest concentration of closed or periodic orbits 
in all panels shown in Figure ^ is near the origin where 
u — v — and orbits are nearly circular. When the spiral 
pattern nears the 2:1 ILR (top panel with ft s = 0.4), the 
periodic orbits become elongated (highly elliptical) and 
so their position in the u, v plane is shifted away from 
the origin. The orbits are elliptical and so have non- 
circular velocity components. On the lower part of this 
panel stars in orbits perpendicular to those supporting 
the arm can be seen. Stars on these orbits have guiding 
radii on the other side of the 2:1 ILR. 

In these panels features are also seen at locations sep- 
arate from the largest concentration at the origin. In 
the second panel from the bottom, structure is seen near 
v ~ 0.05. In this simulation the Sun is just inside the 4:1 
ILR. These orbits correspond to square shaped orbits, 
elongated because of this resonance. To illustrate these 
we show accompanying orbits for these locations in the 
u,v plane. In Figure orbits associated with two of the 
dark regions on this panel are shown. The structure at 
v ~ 0.15 in the panel on figure ^ corresponds to triangu- 
lar shaped orbits influenced by the 3:1 ILR. This orbit is 
also shown in Figure |21 

Nearer to corotation (the lowest panel in^at ft s = 0.7) 
orbits above and below the central family are associated 
with orbits on cither side of the 4:1 resonance (support- 
ing and opposite to the pattern). At higher v, higher to 
Lindblad resonances are encountered closer to the coro- 
tation resonance. Here many resonances influence the 
weighting function and so the model stellar velocity dis- 
tribution. 

4.1. Near the J^:l Inner Lindblad resonance 

In the previous section we discussed two-armed pat- 
terns over a large range of pattern speeds, but at only 
one angular offset, 72 = 45°. We noted in Figure Q 
that at higher pattern speeds more structure appears 
in the weighting function on the u, v plane as higher 
m-Lindblad Resonances are encountered closer to coro- 
tation. In this section we consider in more detail the 
possible role of the 4:1 Inner Lindblad resonance. Be- 
cause square shaped orbits are excited near this reso- 
nance we expect that the velocity distribution could de- 
pend strongly on the orientation of the spiral pattern. 
Therefore we must consider the sensitivity as a function 
of the angular offset 72. 

In Figure weighting on the u, v plane is shown for a 
somewhat stronger two-armed spiral pattern. Parame- 
ters for the spiral pattern are listed in Tabled From top 
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to bottom on these panels, the angular offset between 
the spiral maximum and the Sun 72 = 15, 30, 45, 60 and 
75°. From left to right on this figure the spiral pattern 
speed n s = 0.60,0.625,0.65, and 0.675. In Figure we 
note that a wealth of structure exists in these plots par- 
ticularly for the pattern speeds near and above the 4:1 
ILR at n s « 0.65. 

If the spiral pattern places the Sun near its 4:1 ILR, we 
see that two dominant clumps in the weighting function 
appear near the origin in the u, v plane for low values of 
72 (top panels on rightmost two columns in Figure |3J) . 
In Figure 0] we show the shape of closed orbits (in the 
frame rotating with the pattern) for these two clumps for 
72 = 15° and £l s = 0.675 (also see Table Orbits are 
plotted for (u,v) = (—0.065, —0.043), the inner diamond 
orbit and (u,v) — (0.025,0.015), the outer rectangular 
orbit. The inner orbits would have density peaks at the 
location of two of the density peaks aligned with the 
density perturbation. Consequently, we say they support 
the spiral pattern. The outer orbits have guiding radius 
on the other side of the 4:1 ILR and so are out of phase 
with the inner orbits. 

As pointed out by previous works, two- armed den- 
sity perturbations can excit e square shaped orbits 
ijContopo ulos Sz Grosbol l986f) and a four armed gaseous 
response (Martos ct al. 2005) near the 4:1 ILR. At guid- 
ing radii exterior to the 4:1 ILR, a pure m — 4 poten- 
tial perturbation (with potential components A± 7^ 0, 
A2 = 0), excites orbits that support the spiral pattern, 
Interior to the 4:1 ILR, the closed orbits would be out 
of phase with the density perturbation and so would fail 
to support the spiral pattern. The orientation or phase 
of the orbits changes at the resonance. When the orbits 
support the spiral pattern, they lie on top of the pertur- 
bation in such a way that a self consistent model can be 
created. In other words the density perturbation causes 
orbit perturbations which in turn are consistent with the 
density perturbation. When there is only a strong two- 
armed potential perturbation, as is the case currently 
discussed here, the orbits can still be strongly influenced 
by the 4:1 ILR. In this case the 4:1 ILR is second or- 
der in the epicyclic amplitud e (cx j rather than j 1 ^ 2 ; 
iContopoulos &: Grosbo]|[r986j) . Again the closed orbits 
can be square shaped near the 4:1 ILR, but are oriented 
so that they support the spiral pattern interior to the 
4:1ILR and fail to support the pattern exterior to the 
4:1 ILR. The tendency for orbits to fail to support spi- 
ral str ucture outside the 4:1 ILR prompted iContopoulosI 
(1985) to propose that two-ar med spiral patterns end 
near their 4:1 ILRs (also see IContopoulos fc Grosboll 
ll986UPatsis. k Kaufmarmlll999fl . 

We now compare the structure seen in Figure OH 
with the velocity components of moving groups. In 
our units dominant moving groups in the solar neigh- 
borhood velocity distribution have velocity compo- 
nents (u,v) — (—0.05,-0.10) for the Pleiades mov- 
ing group, (—0.18, —0.09) for the Hyades moving group, 
(—0.04, —0.02) for the Coma Berenices moving group, 
and (0.04, 0.01) for the Sirius/Ursa Major moving group. 
These velocity components are based on the stellar dis- 
tribution measured with Hipparco s observa t ions a nd are 
given wi th respect to the Su n bv iDehnenl 1^998). The 
study bv lFamaev et al.1 l)2004l) grouped the Pleiades and 



Hyades moving groups together and identified the Coma 
Berenices and Sirius/UMA groups as separate structures. 
Note we have divided the velocities by Vb- 

Based on the rightmost two columns in Figure [3] which 
are at pattern speeds £l s = 0.65, 0.675 near the 4:1 IRL, 
it is attractive to associate the Pleiades/Hyades mov- 
ing group with the nearly closed orbits at u ~ —0.05 
and the Coma Berenices moving group with the nearly 
closed orbits near the origin, (u ~ v ~ 0). For 72 ~ 15°, 
the u, v weighting function shows an extension to low u 
for the lower clump, suggesting that low oscillation or- 
bits exist over an elliptical region in the u, v plane that 
could encompass both the Hyades and Pleiades moving 
groups. If this choice of spiral perturbation accounts for 
the Hyades/Pleiades and Coma Berenices moving groups 
then the Hyades/Pleiades moving group stars are in or- 
bits oscillating about the interior diamond shaped closed 
orbit passing through the solar neighborhood (see Figure 
0J and the Coma Berenices group stars are near the outer 
rectangular shaped closed orbit. We note that both fami- 
lies of closed orbits are elongated or not perfectly square. 
This would be expected because the potential perturba- 
tion is elliptical or two-armed. The square shape of these 
closed orbits is a result of the proximity of the 4:1 ILR. 

We have found that the 4:1 ILR provides a promis- 
ing explanation for some of the structure in the solar 
neighborhood velocity distribution. If a spiral density 
wave can produce this structure, then it should be con- 
sistent with the location of spiral arms near the Sun. We 
now consider the location of the spiral arms for a model 
that could account for the Hyades/Pleiades and Coma 
Berenices moving groups. Of the panels shown in Figure 
13 we choose that with fl s = 0.675 and 72 = 15°. The po- 
sition of the two spiral arm density peaks consistent with 
this model are shown as solid lines in Figure 03 If there 
were an additional two arms located in between the dom- 
inant two stellar ones, they would be located where the 
dotted lines are shown on Figure |3J Axes on this figure 
are given in kpc so that the spiral arms in the model can 
be compared with the l ocation of ob served arms (based 
on names referred to by Vallee 2002). 

From Figure |5] we see that our spiral model success- 
fully places the Sun between the Perseus and Sagittar- 
ius/Carina arms. We have adjusted the pitch angle of 
the spiral arms (set a) so that the separation between 
the Perseus arm and the Sagitt arius/Carina arm is con- 
sistent with that suggested bv IValleel l)2002fl . However, 
our model requires that the Sun is closer to the Perseus 
arm than the likely location of the Sagitta r ius/ C arina 
arm, in contradiction to the model bv IValleel l|2002fl . We 
note that the Sagittari us/Carina arm is n o t as s trong 
in stars as it is in gas Drimmcl & Spcrgcl (2001), and 
the dust map of Drimmcl & Spcrgcl ( 2001) shows a kink 
or bend in the Sagittarius/Carina near the Sun. This 
bend may be associated with the ch ange in the o rbit 
structure at the 4:1ILR. The model of I Vailed 1)20021) as- 
sume purely a logarithmic shape for the spiral arms, 
whereas the orbit structure near the 4:1 probably causes 
deviations in the location of the arms. Fourier decom- 
position of nearby galaxy images has shown that even 
when the galaxy is grand design, Fo urier components 
with to > 2 may still be significant (e.g.. lElmeg reen et all 
H992tlSeigar fc Jameslll998HGrosbol et al.ll2004|) and the 
simple model explored here only considers a dominant 
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m = 2 Fourier component. A more detailed model for 
the arms (based on the observations) and a self consistent 
model for the spiral arms, the gas response and popula- 
tion of orbits would be required to determine if the arm 
locations were consistent with the proposed dynamical 
model. 

We now discuss sensitivity of our proposed spiral wave 
model to the parameters used to define it. We have car- 
ried out comparison simulations with different spiral arm 
pitch angles. We find that the weighting function con- 
tours (and inferred structure in the velocity distribution) 
are not strongly sensitive to the value of ai. For ai 
between 5-7 , consistent with th e ran ge of p itch angles 
discussed bv lDrimmel fc Spergelll200HlValleeil2002l little 
difference is seen the structure of the weighting function. 
However, from Figure [3] it is clear that the model veloc- 
ity distribution (and orbital structure) is very sensitive 
to the orientation angle of the spiral arms (72) with re- 
spect to the Sun. Changes in 72 of as small as 10° can 
cause significant changes in the locations of closed orbits 
in the solar neighborhood. 

We have also varied the strength of the two-armed per- 
turbation, A2. For \A2 1 < 0.005 the orbits are more 
nearly circular and would not be able to be consistent 
with two dominant spiral arms (the orbits are not suffi- 
ciently elliptical away from the 2:1 ILR). For \ A 2 \ < 0.005 
the width of the 4:1 ILR is smaller (the resonance is sec- 
ond order the epicyclic amplitude) and the two families 
of orbits near the 4:1 ILR are closer together, too close 
to account for the Hyades/Pleiades and Coma Berenices 
moving group separations. For IA2I > 0.01 the sepa- 
ration between the two orbit families was larger than 
that of the two moving groups. We had best results 
finding a model u, v distribution similar to that ob- 
served in the solar neighborhood for A2 ~ 0.008 ± 0.002. 
We note that this is about twice a s high as that esti- 
mated from the K-ba nd model bv IDrimmel fc Spergell 
l)2001() . Specifically IDrimmel fc Spergell (|2001f) esti- 
mated (T, max - S min )/E mm = 0.32 which is smaller 
tha n the contrast seen in other similar near by galax- 
ies llSeigar fc Jameslll998t iG rosbol et al.ll2004l) . Conse- 
auentlv lDrimmel fc Spergell l)2001| ) commented that they 
suspected the true value for the density contrast (arm vs 
interarm) could be higher than this number. 

In Figure we directly compare our model (see Table 
^for parameters) it, v plane panel showing our weighting 
function, to the observed stellar velocity distribution in 
the solar neighborhood. The middle panel on this figure 
is the stellar v elocity distribution of th e F and G dwarf 
stars listed bv iNordstrom et all l|2004ft . The righ t hand 
panel of this figure is from Figure 3 of lDehnenl lfl99l 
which shows all s tars in his Hipp arcos sample. This panel 
from Figure 3 bylDehner] l)1998f) resemble s Figure 20 by 
INordstrom et all ( 2004) and Fi gure 7 from Famaev et al. 
(2004). In our panel the origin corresponds to a circu- 
lar orbit in the absence of perturbations by spiral arms, 
however the stellar velocity distribution is plotted with 
respect to the solar velocity. If our proposed dynamical 
model is correct it suggests that the Sun is moving with 
u,v > 0. 

In the Figure bv lDehnenl l)1998l) the clump nearest the 
origin is the Coma Berenices moving group. The clump 
below the Coma Berenices group is the Pleiades mov- 
ing group. The clump at positive velocities is the Sir- 



ius/UMA moving group and the clump at v — — 40km/s 
is the Hyades moving group. We note that our model 
predicts that the Hyades and Pleiades moving groups 
are kinem atically related. Th i s is consistent with the 
studies bv lFamaev et alJ l|2004D : INordstrom et all ( 2004) 
which found no strong separation between them. Our 
model is consistent with the velocity separation between 
the Coma Berenices and Hyades/Pleiades moving groups 
and predicts a slight tilt (larger v for more negative u) to 
the Hyades/Pleiades moving group contours that is seen 
in the observed velocity distribution. 

Based on the comparison between our model and the 
observed velocity distribution, it is tempting to asso- 
ciate the Sirius/UMA moving group with closed orbits 
at higher v that would correspond to higher m ILRs. 
However since these would be affected by higher m res- 
onances, the structure of the orbits is likely to be very 
strongly dependent on the location of the Sun with re- 
spect to the spiral structure and on the assumed struc- 
ture of the spiral arms. Consequently we don't feel we 
can attach any significant to structure in our model at 
the location of the Sirius/UMA moving group. If higher 
order resonances affect stars with guiding radius just out- 
side the Sun, then it's possible t hat the outer Milky Way 
disk is flocculent ( Quillen 2002) . It is also possible that 
an additional spiral density wave could affect the outer 
part of the galaxy, adding an additional complication 
that we have not explored here. 

The 4:1 Lindblad resonance provides an explanation 
for some of the structure present in the stellar velocity 
distribution in the solar neighborhood. However, this 
explanation succeeds only over a very narrow range of 
pattern speeds. For the small range of reasonable pertur- 
bation strengths, based on the structure seen in FigureEl 
we estimate that the spiral pattern speed must be within 
fl s = 0.66 ± 0.03. This corresponds to a pattern speed of 
ft s = 18.1±0.8 km s _1 kpc _1 (after restoring the physical 
units). This p attern speed is consistent with many previ- 
ous estimates (lShavivl2003l) and the k inematic models by 
iMartos et all <|2005fk Bissantz et alJ l|2003D . Our model 
for the veloc i ty dis tribution is quite similar to that of 
IMartos et alJ l)2005f) who showed that a dominant two- 
armed stellar perturbation could cause four arms in the 
gas distribution. 

5. SUMMARY AND DISCUSSION 

In this paper we have considered the effect of spiral 
density waves on structure in the stellar velocity distri- 
bution in the solar neighborhood. We find that Lindblad 
resonances with spiral density waves can cause structure 
in the velocity distribution. Lindblad resonances can ex- 
cite large epicyclic perturbations, allowing stars from dis- 
tant locations to reach the solar neighborhood. In the so- 
lar neighborhood, we find that the location of a resonance 
is primarily set by the tangential velocity component, v. 
Because spiral perturbations are not strong, their widths 
are narrow and they only excite large epicylic oscillations 
over a narrow range of guiding or mean radii. However 
because these epicylic oscillations can cause stars to cross 
into the solar neighborhood, they can still cause struc- 
ture in the solar neighborhood velocity distribution. 

Because Lindblad resonances can cause significant 
changes in the structure of orbits, they provide a promis- 
ing explanations for structure seen in the solar neighbor- 
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hood velocity distribution, or moving groups and super- 
clusters. To explore this possibility we have searched for 
orbits that are near closed or periodic orbits in the frame 
moving with the spiral pattern. We have constructed a 
weighting function that estimates the distance in phase 
space of an orbit from a closed or periodic orbit. When 
a moderately strong two-armed perturbation is present 
with pattern speed at 0.66 ± 0.03 times the angular ro- 
tation rate of the Sun, placing the Sun near the 4:1 
Lindblad resonance, two regions near the origin in the 
u,v plane exhibit nearly closed orbits. These islands are 
likely to be populated by stars spanning a range of ages, 
as are moving groups. 

For a two-armed model with pattern speed placing the 
Sun near the 4:1 Lindblad resonance, and an angular 
offset of 72 = 15°, the structure of our weighting func- 
tion resembles the observed velocity distribution. Nearly 
closed orbits are located near existing moving groups. 
The region at at negative v in the u, v plane is elon- 
gated toward negative u and so could encompass both 
the Hyades and Pleiades moving groups. The region at 
positive v, can be associated with the Coma Berenices 
moving group. This model provides a possible expla- 
nation for some of the structure observed in the solar 
neighborhood's velocity distribution. The model suc- 
ceeds only over a very narrow range of pattern speeds; 
with fl s = 18.1 ± 0.8 km s^ 1 kpc -1 providing a tight 
constraint on the angular rotation rate of the spiral pat- 
tern. Our model is m ost similar to and co nsistent with 
the model proposed bv lMartos et alJ l)2005l) in which two 
dominant stellar arms excite a four-armed gaseous re- 
sponse. 

We have shown here that resonances with spiral pat- 
terns provide a promising way to explain clumps in the 
solar neighborhood velocity distribution. In this paper 
we have only considered the response of stars to spi- 
ral density perturbations but did not consider the re- 
sponse of the gas. In our model we have not constructed 
a stellar orbit distribution consistent with the assumed 
density perturbation. Future work should strive to cre- 
ate models that are more self-consistent. Here we have 
not explored the sensitivity of the ste llar distr i butio n 
with stellar birth site or age. However IPehnenl {l998); 
iFamaev et all l)2004|) found that structure in the velocity 
distribution depends on the type of star or stellar popu- 
lation. Here we used a crude weighting function to find 
nearly periodic orbits in the u, v plane, however we have 
not produced model stellar density distributions. Future 
work could consider the birth of stars in a disk supporting 
a spiral pattern and explore ways to predict the location 
and number of stars as a function of birth site and age. 
Future work could also consider the sensitivity of the 
present day velocity distribution to the way that spiral 
structure evolves. For example it may be possible to dif- 
fer entiate between pure am plitude growth (as considered 
by iDe Simone et all l2004|) from shearing density wave 
models for which both amplitude an d wavevector v ary 
simultaneously (e.g., as considered bv lFuchsll2001albfl . 

There are some interesting consequences of our model. 
The 4:1 resonance is likely to cause deviations from a 
pure logarithmic spiral pattern near the Sun. Better 
models for the Galaxy could include arms that deviate 
from logarithmic spirals. Because the Sun may be lo- 
cated near a Lindblad resonance causing large epicyclic 



amplitudes in the stars, measurements of the velocity of 
the local standard of rest (VLSR) and Oort's constants 
are likely to be biased. By understanding and correct- 
ing for these biases, measurements of these astronomical 
quantities may be improved. In this paper we have not 
proposed a dynamical explanation for the Sirius/UMA 
moving group. This group could be related to a higher 
m Lindblad resonance, or perhaps there is another spi- 
ral density wave moving at a slower pattern speed at 
larger Galactocentric radius. Here we have suggested 
that the Hyades/Pleiades moving groups are on orbits 
with mean radii within the Sun's Galactocentric radius, 
and the Coma Berenices and UMA/Sirius moving groups 
are on orbits with mean or guiding radii outside Rq. It 
would be interesting to see if these stellar populations 
have different age and metallicity distributions. Here we 
have not discussed the orbits of stars at velocities in be- 
tween the moving groups. Because of the 4:1 ILR, these 
stars should have large epicylic amplitudes and so could 
be part of a different stellar population than the moving 
groups. 

Here we have not considered the role of more than one 
spiral density wave or the role of the Galactic Bar. Stel- 
lar orbits that are not affected by a Lindblad resonance 
from an additional perturbation, such as from another 
spiral density wave or the Galactic Bar, would be only 
weakly affected by the additional perturbation. Away 
from a resonance, the orbital kinematics can be treated 
wit h low order perturbation t heory (for example as done 
bv iBinnev & Tremaiiiel 11987ft . Only stars (at v ~ 50 
km/s) associated with the Hercules stream that are in- 
fluenced by the Bar's 2:1 Lindblad resonance, should be 
strongly affected by the Galactic Bar. Stars at lower 
values of \u\, \ v\, should be distant and so unaffected by 
Lindblad resonances with the Galactic Bar. We expect 
that the closed orbits considered in this paper weakly os- 
cillate or are weakly perturbed at the frequency of the 
Galactic Bar. These oscillations should cause small vari- 
ations in the location of the periodic orbits on the u, v 
plane as seen from the Solar neighborhood. However, as 
long as we consider orbits that are not associated with 
resonances with the Galactic bar, then strong features 
associated with the Bar would not be seen in the veloc- 
ity distribution. The same situation is likely if there is 
an additional spiral density wave (at a different pattern 
speed) present i n the solar neighborh ood (as suggested 
most recently bv lNaoz fc Shavivll2005ft . We only expect 
strong structures in the velocity distribution that are as- 
sociated with resonances from one of the spiral density 
waves present in the Solar neighborhood. However, the 
presence of more than one perturbation can influence the 
stellar dynamics. For example, if the Solar neighborhood 
is affected by more than one spiral density wave, then 
stars at veloci ties between moving groups might be on 
chaotic orbits l|Quillenl l2003). a factor which could cause 
a relatively large increase in their velocity dispersion with 
time. 
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n=2 7=45 A=-0.0050 a=7.0 Q =0.40 



n=2 7=45 A=-0.0050 a=7.0 0=0.50 



n=2 7=45 i=-0.0050 .7 = 7 0=0.60 




77=2 7=45 A=-0.0050 r=7.0 5 =0.70 




Fig. 1. — The (u, v) plane shown for two-armed spiral models at different pattern speeds. The u, v axes are shown in units of Vb, the 
velocity of a star in a circular orbit at the Sun's Galactocentric radius. The gray scale and contours show our weighting function W which 
gives a measure of the extent of epicylic motion about a closed or periodic orbit. From top panel to bottom panel we show the effect of 
spiral structure with pattern speed Q s = 0.4,0.5,0.6 and 0.7 respectively (in units of Vo/Ro)- The two armed perturbations have angular 
offset 72 = 45°, wavevector with a = 7, and potential perturbation strength A2 = 0.005 (see Table0. Dark regions correspond to orbits 
which are nearly periodic. The top panel is nearest the 2:1 ILR and so closed or periodic orbits (in the frame of the pattern) are highly 
elliptical. This causes the periodic orbits to be moderately distant from the origin where u = v = and where circular orbits are located 
in the absence of spiral perturbations. The bottom panel is near the outside the 4:1 Lindblad resonance. The structure at low v in the 
bottom panel is caused by the 4:1 Lindblad resonance. At larger v structure from higher m Lindblad resonances is encountered closer to 
corotation. 
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Fig. 2. — Closed orbits in the frame moving with the spiral structure are shown for the two-armed Q s = 0.60 model (see Table IT| for 
parameters). The weighting function for this model is shown as the second panel from the bottom in Figure Orbits are shown for two 
dark regions shown on this panel. The triangular closed orbit corresponds to position on the u,v plane with u = 0.00, v = —0.15. The 
square orbit corresponds to u = —0.20, v = 0.05. 
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Fig. 3. — The (u, v) plane shown for stronger two-armed spiral models at pattern speeds placing the Sun near the 4:1 Inner Lindblad 
resonance. Panels are similar to those shown in Figure From the leftmost column to the rightmost the pattern speed is Q a = 
0.60, 0.625, 0.65 and 0.675 respectively. In each column from top to bottom the angle of the spiral pattern with respect to the Sun is 
7s = 15,30,45,60 and 75°, respectively. Parameters for these simulations are also listed in Table HI 
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Fig. 4. — Closed orbits in the frame moving with the spiral structure are shown for the two-armed Q s = 0.675 model (see Table U for 
parameters). The weighting function for this model is shown as the top panel in the third column from left in Figure [3] Orbits are shown 
for two dark regions shown on this panel at the location of the Sun. The inner square closed orbit corresponds to position on the u, v plane 
with u = —0.065, v = —0.043. This orbit supports (lies on top of) the dominant two-armed perturbation. In other words, two of the orbit 
peaks are on top of two of the arms. The outer rounder square orbit corresponds to u = 0.025, v = 0.015 and is out of phase with the 
two-armed density perturbation. We also show an inner orbit supporting the two-arms inside the Galactocentric radius of the Sun. 
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Fig. 5. — Location of arms with respect to the Sun consistent with the distributions shown in Figure El with 72 = 15° (top panels). The 
two stellar arms are shown as solid lines. If two additional arms were present in between the two strong stellar ones, they would be located 
approximately at the dashed lines. The position of the Sun is shown as a small circle. The location of the 2:1 ILR, 4:1 ILR and corotation 
resonance are shown as large dotted circles. The angle between the Sun and nearest strong arm, 721 is that between the two vectors. We 
have labeled the arms according to their common names. In this figure Galactic rotation is clockwise. See Table HI for descriptions of the 
parameters. 
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Fig. 6. — Comparison between our model for the u,v plane (leftmost panel) and the observed stellar velocity distribution. The left 
hand panel (similar to Figure ^ has parameter s listed in Table HI The middle panel is a 2D histogram of the u, v velocities of the F and 
G dwarf stars listed bv N ordstrom et alJ 120041) in units of Vb- The u, V velocity components listed in Table 1 by [Nordstrom ct al. 120041) 
are heliocentric. To compare these velocities to those shown in the lefthand panel we have assumed a VLSR velocity of Uo ~ Vb ~ 8km/s. 
The right hand panel is from|Dchncn 1 199£j) Figure 3 for all stars in his Hipparcos subsample. In the right panel, the clump nearest the 
origin is the Coma Berenices moving group. The clump below the Coma Berenices group is the Pleiades moving group. The clump at 
positive velocities is the Sirius/UMA moving group and the clump at v = — 40km/s is the Hyades moving group. We note that our model 
predicts that the Hyades and Pleiades moving groups are kinematically related. This is consistent with the study by fFamacy ct aLl l2004|) 
which found no strong separation between them. Our model is consistent with the velocity separation between the Coma Berenices and 
Hyades/Pleiades moving groups and predicts a slight tilt (larger v for more negative u) to the Hyades/Pleiades moving groups contours 
that is seen in the observed velocity distribution. In our model, the Hyades/Pleiades stars are on diamond shaped orbits supporting the 
spiral pattern and those in the Coma Berenices group are on orbits that are out of phase by 45°. The Hyades/Pleiades stars have mean 
radii that are within the Sun's Galactocentric radius. Those at larger v have mean radii outside Rq. 
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TABLE 1 

Parameters describing spiral patterns 



Figures 


A 2 


72 


OL2 


n s 


1 
2 
3 
1 


5l6l 


-0.005 
-0.005 
-0.008 
-0.008 


45° 
45° 

15,30,45,60,75° 
15° 


-7 
-7 
-6 
-6 


0.4,0.5,0.6,0.7 
0.6 

0.60,0.625,0.65,0.675 
0.675 



Note. — The parameters of the spiral pattern corresponding 
to simulations shown in the Figures. The angular offset 72 in 
degrees are the angle (measured at the Galactic center) between 
the location of the Sun and the nearest density maximum of the 
spiral potential perturbation at a radius of Ho (see Figure^- The 
perturbation strength A2 is given in units of Vq the velocity of 
a star in a circular orbit at i?o- The pattern speed, Q s , is given 
in units of fig = Vo/Rq- The parameter 0L2 sets the pitch angle 
of the spiral arms. For 02 < the arms are trailing when the 
rotation is clockwise. 



